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Abstract 

In acoustic, ultrasonic or electromagnetic propagation, crossed me- 
dia are often modelled by linear filters with complex gains in accor- 
dance with the Beer-Lambert law. This paper addresses the problem of 
propagation in media where polarization has to be taken into account. 
Because waves are now bi-dimensional, an unique filter is not sufficient 
to represent the effects of the medium. We propose a model which 
uses four linear invariant filters, which allows to take into account ex- 
changes between components of the field. We call it bi-filter because 
it has two inputs and two outputs. Such a circuit can be fitted to light 
devices like polarizers, rotators and compensators and to propagation 
in free space. We give a generalization of the Beer-Lambert law which 
can be reduced to the usual one in some cases and which justifies the 
proposed model for propagation of electromagnic beams in continuous 
media. 

keywords: linear filtering, polarization, Beer-Lambert law, random 
processes. 

1 Introduction 

Propagation of acoustic or electromagnetic beams is often studied through 
linear differential equations with coefficients which depend on the medium 
characteristics. Because of the difficulty to estimate parameters and to give 
well-fitted equations, the medium is often taken to be a linear invariant 
filter (LIF) with given spectral gains. For instance, propagation of ultra- 
sonics through biological tissues is taken to be a filter with complex gain 
exp[— (W 3 ] . In scanning, values of parameters allow diagnosis [3], [10], [7]. 
The same model is used for losses and dispersion of radiowaves in coaxial 
cables, free space or fiber optics [5], [S]. 

We consider beams reduced to a trajectory in some medium. Beams 
measurements are very often reduced to the values of quantities available 
at two points, one at the origin and one at the end point (for instance 
the amplitude or the power). Even if the beam is received in a antenna 



with a given area, the information is often added in a coaxial cable or in a 
wave guide, and reduced to a finite number of complex quantities. Moreover, 
comparisons about values at the transmitter and at the receiver give insights 
about the medium, and it is the aim of most devices used in practice (and 
this is true for waves in any frequency band, optics, radio, acoustics or 
ultrasonics). Actually, the main difference between sonics or ultrasonics 
waves and radio or optical waves is the polarization of the last one. Then 
they are defined by vectors instead of scalars. This paper highlights this 
difference. 



We define a one-dimensional beam by an electric field E =<^E (t) , t € R j 
at each point of a trajectory, t stands for the time, z is the coordinate on 
the trajectory and E (t) is a vector orthogonal to the trajectory. In nu- 
merous situations, it is sufficient to consider the amplitude of E^ (t) , which 
defines the power and the spectral content of the wave. Then, it is suffi- 
cient to model the medium as a linear filtering, as in the case of acoustic 
wave propagation. It is no longer the case when polarization phenomena 
have to be taken into account. We have to show the evolution of two com- 
ponents E z x it) , Ey (t) of the electric field which defines the wave. Because 
both components can be linked, the behavior of one component at the time 
t can depend on the other component. In a linear model, this means that 
each component is the addition of linear filtering of itself and of the other 
component. Consequently, the medium has to be defined by a family of 
four filters. Nevertheless, we will talk about a bi-filter because the opera- 
tion defines a linear way between two couples of processes, two inputs and 
two outputs. We place ourselves in a stationary frame where the processes 
are stationary (in the wide sense) with stationary correlations. The filters 
which model the medium will be time invariant with the acronym LIF for 
Linear (Time) Invariant Filters and will be defined by complex gains rather 
than impulse responses. They generalize the notion of "scattering matrix" 
used for monochromatic waves where each coordinate of the scattered wave 
is a linear combination of transmitted coordinates. Though power spectra 
of electromagnetic waves are the most often assumed in limited bands, the 
definitions and computations will be done for general spectra possibly with 
infinite support. 

Section 1 of this paper addresses the definition and properties of bi- 
filters. Sections 2 and 3 apply results of section 1 to one-dimensional beams. 
Section 4 studies a generalization of the Beer-Lambert law which shows why 
a continuous medium can be modelled by a bi-filter. Appendix gives proofs 
of formulae. 
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2 Bi-filtering 

1) We consider two (real or complex) stationary processes 

Xi = {Xi (t) , t € M} , Yi = {Yy (t),t€ M} 

with spectral density sx 1 , «Yi depending on the frequency lo/2tt, but we will 
omit the variable u when the result is not ambiguous. Moreover we assume 
a stationary correlation between them which defines a cross-spectrum sxiY a 
such that P, [12] 

E [X x (t) X* (f - r)] = SXl (u) e^du 
E [y (i) Y* (t - r)] = fZo s Yl (w) e^dw (1) 
E [X x (t) (t - r)] = SXlYl H e^dw 

where E[..] stands for the mathematical expectation (or ensemble mean) and 
the superscript * stands for the complex conjugate. 

Now, we define the processes X 2 ={X 2 (*),*€ M} ,Y 2 =-jT 2 (t),iG K} by 

X 2 (i) = X x * /i n (t) + y * /i 12 (t) 

y 2 (t) = x x * h 21 (t) + y * h 22 (t) 1 J 

where the hjk can be considered as impulse responses of four linear invariant 
filters Hjk and (. * .) stands for the convolution product. For example we 
have 

/oo 
X\ (u) h 2 \ (t — u) du 
-oo 

We will say that (Xi,Yi) and (X 2 ,Y 2 ) are the input and the output of 
the bi-filter T-L = {Tijk,j, k = 1, 2} . We know that it is more convenient to 
use spectral gains Hjk (uj) rather than impulse responses. When impulse 
responses are sufficiently regular we have together 

(*) = h I-oo H 3k M e^du 
H jk (u) = J^ 00 h jk (t)e-^ t dt {6) 

but the hjk are not always ordinary functions. The following writing of (2) 
is more general because the impulse responses do not appear 

X 2 (i) = ftn[Xi](i) + ^i2[Yi](t) 

y 2 (t) = h 2 i [Xi] (t) + % 22 [Yi] (t) y > 

It is proved in appendix 1 that the bi-dimensional process (X 2 ,Y 2 ) is sta- 
tionary with spectral characteristics perfectly defined by formulae (34) which 
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Figure 1: 



are simplified in 

S X 2 Y 2 = #12#22 s n + HiiHl x s Xl + Hi 2 H2 1 Sy 1 X 1 + #ll#22 s YiXi 

s X2 = \H 12 \ 2 syi + \H U \ 2 s Xl + 2TZ [H^H^syiXi] 
sy 2 = \H 2l \ 2 s Xl + \H 22 \ 2 s Yl + 2K [H 22 H^s YlXl } 
< s Xl x 2 = H$ 2 s YiXi + H* lSXl (5) 
s Xl y 2 = H^ 2 s YiXi + H^s Xl 
«Fii 2 = H^SYiXi + H{ 2 s Yl 
SYiy 2 = H 2l s YlXl + H 2 * 2 s Yl 

where 1Z [..] stands for the real part, s a and s a f, are for spectra or cross- 
spectra, Hjk for filters complex gains, and all terms depend on uj. The 
figure 1 gives the scheme of the bi-filter. When used in the matrix form 
H=[Hij] , the bi-filter is a "scattering matrix" for each value of uj. 

3 Application to waves 
3.1 General formulae 

1) We consider a beam which crosses a medium along the axis Oz of the 
orthogonal trihedron Oxyz. The beam is defined by its electrical field 
~Ej (t) = (E* (t) , Ey (i)) at time t and distance z, where components (which 
are orthogonal to Oz) are taken on the axes Ox and Oy. We assume 
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a) that the processes (at O) E° = {E® (t) , t G R} and E° = {E® (t) , t G 1 
are stationary and with stationary correlation. The spectra and the cross- 
spectrum are 



si M 



- 4 M , (w) 



i?J^ 2 (w) such that, whatever < u < z 



(6) 



EZ(t) = wg\E«\ W + ^f [e«] (t) 

This writing highlights the dependence of E z (t) and E z (t) on the whole 
processes E" and E^ and not only on the r.v. EJ£ (t) and Ey (t) . Formulae 
(5) allow to write the different spectra and cross-spectra of the beam as 
(with Hjk = Hjj; to alleviate formulae) 



I tj 1 2 „0 i I u |2 



H 12 H2 1 s yx + HnH^s^y 



s% = \H l2 \ 2 s° y + \H n \ 2 4 + 2^ [^12^1*1^] 

S + l^2 2 | 2 4 + 2^[ J ff 22 i/ 2 * lS y 



2/ 

-02 

'rrrr 

02 

-01 



\H 2 i\'i 

H^sfy + H. 
H* n s° vx + H* 12 s° v 







22 s y 



} 



b) and that the medium between the points z\ and z 2 can be charac- 
terized by four LIF (linear invariant filters) H 2 ^? 2 with the complex gains 



(7) 



where the seven equalities are respectively for the couples (from the top to 
the bottom) 



(E z E z ) (E z E z ) ( E z E z ) ( E° E z ) (E° E z ) (E° E z ) (E° E z ) 
y-'xi ) ' K^xi ^x) i V y y ) > V xi ^x) > y-'xi ^y) > V y x) > V y y) ■ 

The power P z of the wave at z is defined from the components by 



E 



\E Z (t)\ 2 +\E Z (t)\' 



\_ s x s y\ 



(lo) dco. 



(8) 



This definition is justified (at a multiplicative constant in accordance with 
some physical system of units) because devices for intensity measurements 
are not sensitive to the direction of the electric field. Also, the definition is 
independent of the basis used (see below). 

2) In the proposed model, the filters complex gains are proper to the 
coordinate system Oxy. For the system Ox'y' deduced by a rotation of 
angle 6, the new system of complex gains K z jf 2 is defined by (we give the 
system for = K jk ) 







" H u - 




= 


H12 


K 2 i 


H21 


K22 




H22 
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p cos e 

-P sin 9 



Psin6> 
Pcos(9 



cost/ 
- sin( 



sm( 
cos I 



(9) 



is orthogonal and then 1 = 0* (the transpose of 0). Also, the only 
set of parameters which are invariant by rotation verify (see appendix 2) 



H 



ii 



22- 



H 



12 



-H- 



21- 



(10) 



This set of Hjj~ defines the proper subspace of the eigenvalue 1. This prop- 
erty is very important and will often be used. Moreover it is easy to verify 
that 

Hn + #22 = Kn + K 2 2 ,j-q 

H11H22 — H12H21 = K\\K 2 2 — Ki 2 K 2 \. 

This property of invariance will be used in section 4.2. 

3) The sum of bi-filters /C = % + H' is naturally defined by Kjk = Hjk + 
H'- k whatever j, k. It is the same as the usual sum for linear filters (filters in 
parallel of the circuit theory) . The product of bi-filters /C = TixH' is defined 
in the same way as the filters in series of the circuit theory 

(X 2 ,Y 2 ) = /C[Xi,Yi] =%'{H [Xi,Yi]} 

The scheme of this operation is given in figure 2 with the notations used 
in the next sections. Equivalently, we have (because the complex gains of 
LIF in parallel and LIF in series are the sum and the product of individual 
complex gains) 

Ku = HuHL + H 2 \H{ 



K 12 
K21 
K22 



11 -r "21-"i2 

H\2H'u + H 2 2H' 12 
Hi\H' 21 + H2\H' 22 
H 22 H 22 . 



Hi 2 H' 21 



This operation is associative but not commutative (which is a huge difference 
with ordinary LIF). The filter X such that 

In = -^22 = 1, 112 = hi = 

is the unit: IxH = "HxX = %. 

Sums of bi-filters have to be used when the beam is split so that the 
results cross different media. Each of them corresponds to a bi-filter and 
outputs are added (similar to a fringes pattern). The product is for a beam 
which crosses two successive media, or two successive thickness of the same 
medium, each of them being defined as a bi-filter. 



3.2 Polarized beam 

A "polarized wave" at z is defined by a direction ip (with respect to the axis 
Ox), and some stationary process A Z ={A Z (t) ,t G M} 



E z At) 

E 'v (*) 



A z (t) cos ip 
A z (t) sin-0. 



(12) 
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This definition is consistent because, in the basis Ox'y' such as (Ox, Ox') 
tp' we have 

E z , {t) = A z (t) cos^-V'O 
E z .{t) = A Z (t)sin(^-V')- 



In this paper the "polarized wave" is for the "linear polarized wave" used 
for deterministic beams. Because the sum of two polarized waves (A 2 ,^) 
and (B z , (f>) is generally not polarized, the set of polarized waves has not an 
interesting algebraic structure. But, such waves can be treated separately 
when transformations are linear. It is worth-noting that the wave 

T[E Z X ],T[E Z ] 

where T is any LIF, is polarized in the same direction that (E^E 2 ). 
From (11) , a necessary condition for a polarized wave is 



pIpI 



{pIvY 



(13) 



where 



P. = E \E z (t)\ 



E 



E z (t)\ 2 ],p z xy = E[E z (t)E z *(t)] 



(and not \p% y \ except for real processes). Conversely, if (13) is verified, p 
is real (because p z x p z y > 0), which leads to 



xy 



E 



\E*(t)-\E*(t)f 



Ipl - \e, 



where e = ±1 and A £ C (e = 1 when p xy > 0). Consequently, we have 
E z (t) = XE Z (t) with A = p x / p y (quantity which is real and independent 

of t) and (13) becomes a sufficient condition for polarization. From (12) the 
condition (13) is true whatever the system Oxy. Reciprocally, when verified 
in one particular system, it is verified in the others. Also, the equality (13) 
is true for correlations and then for (regular) spectra 



s z s z 



(Sly) (14) 

but this last equality is not sufficient for polarization. As an example 

f E z {t) =A z {t)cosij} 
\ E z (i) = B z (t)sm4> 

where B z is the output of a LIF with input A 2 and complex gain e (uj) 
taking only the values ±1. The condition (14) is fulfilled but the wave is not 
generally polarized and can be split in the sum of two polarized waves in 
the directions ip and —ip. 



7 



3.3 Unpolarized beam 

1) We define an "unpolarized wave" E^ at z by the condition (whatever the 
systems Oxy, Ox'y') 



q z — e z — n 
°xy °x'y' u- 

Equivalently (see appendix 3) 

q2 — qZ — qZ — „ Z 

By- (15) 
s xy — s x'y' ~ U ' 

The term "unpolarized wave" is used in this paper rather than the term 
"circular polarized wave" which we encounter for deterministic waves. The 
subset of uncorrelated unpolarized waves is a group for the addition. But 
the sum of two correlated unpolarized waves can be not unpolarized. For 
example, (E x ,E y ) added to (E x ,—~E y ) is polarized along Ox. 

In |15j . pp. 350, and in the optics communauty, the unpolarized light 
(natural light) is defined by the equality p z xy = 0, whatever Oxy, which 
implies p z x = p y . This condition is much more weak than the condition 
(15) . In the strong definition, an unpolarized wave remains unpolarized after 
crossing of a compensator. It is not true when using the weak definition. 

2) i? is an unpolarized wave when its components and E^J are un- 
correlated whatever the system Oxy (strong definition) . It is the case if and 
only if the components are uncorrelated for a given system in which the 
spectra are identical (s xy = 0, s x = s y ). In this circumstance, (7) is reduced 
to 



s xy — {H12H22 + H11H21) s x 



ll H 12 \ 2 + \H 11 ^ -° 



s y = ( l-^2i| 2 + \H22\ 2 ) s? 
Consequently the wave E^ remains unpolarized whatever s x if and only if 
H12H22 + H\\H 21 = 

|-ffl2| 2 + |-ffll| 2 = I #21 1 2 + |#22| 2 • 

Equivalently it exists a real function P (oS) such that 

H l2 {oj) = -H* 2l {u)e iP ^. K ' 



3.4 Partially polarized beam 

Any wave which is not polarized or not unpolarized is a partially polarized 
wave. This is a definition which results in different classes following the 
definition given to unpolarization (weak or strong sense). In the class of 
beams defined in appendix 5, a partially polarized beam corresponds to a 
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probability law different of the uniform on (0, 2tt) (case of an unpolarized 
wave) and not degenerate (polarized wave when degenerate). 

The Stokes decomposition theorem is still a studied problem [19]. It 
states that a partially polarized beam is the sum of a polarized beam and of 
an unpolarized beam. For stationary processes with any spectra, solutions of 
the problem can be found according to the set where solutions are searched, 
and according to definition of unpolarization |8]. It is possible to construct 
waves which are polarized in more than one direction. Provided that the 
operations are linear, we can split the beam in a convenient number of 
polarized ones and we can study them separately. 



3.5 Examples of bi- filters 

We characterize the following elementary bi-filters by the matrix H of the 
Hjk (oj) in Oxy or the matrix K of the Kjk (w) in Ox'y': 

_ Hu H12 

The elementary operations which follow are often assumed independent of 
the frequency. Actually this property is generally approached in the (lim- 
ited) frequency band of an experiment. However the matrix H which defines 
the bi-filter is a function of the frequency uj/2it in most cases. 

1) A compensator (9 x ,9 y ) between O and z induces different delays 9 X 
and 9 y for the components: E z x (t) = E%(t- 9 X ) and E* = E° (t - 9 y ) [TH]. 
The equivalent bi-filter is defined by 

_ r e -k>»« 

In the basis Ox'y' defined by ifi =(Ox, Ox'), we have 

K _ e~ tul9x costp e~ luiey sintp ' 



2) The polarizer along the axis Ox' suppresses the orthogonal component 
along Oy'. This defines the bi-filter (in the basis Oxy with ip = (Ox, Ox')) 



H 



cos 2 ip sin ip cos tp 
sin tp cos ip sin 2 ip 



3) The compensator (9 x ,9 y ) followed by a polarizer in the direction Ox' 
defines the bi-filter 



H 



-iuid x 



COS tjj 
sin ip cos if) 



-iu>9 



y sin t/j cos ip 



e -iw9 y gin 2 ^ 



9 



The power P z at z is given by, from (7) and (8) 



pz = IZo K + Pfi + m (74,)] m du 

a = cos 2 ipi P = sin 2 tp 

is a function of #y — X and ^. A good choice of these parameters allows 
to measure particular values of P z . If the parameters do not depend on the 
frequency, this allows to estimate the quantities 



E 

E 
E 



\E° x (t) 
\E° y (t) 



[El it) Ep (t)] = s% (u) 
which define the Stokes parameters and others 115) . 



3.6 Rotator 

In appendix 2, we prove that any polarized beam in the direction ip (with 
respect to Oxy) at O will be polarized in the direction if) + 9 at z (9 ^ 
77 mod 7r) if and only if 

H21 cos 9 = H±i sin 9 

H l2 = -H 21 (17) 
H\\ = H22 

The angle of polarization is changed by the quantity 9. Moreover, the am- 
plitudes at O and z are linked through a LIF with complex gain H\\j cos 9 

A* = J_W U [A°] . (18) 
cos 9 

Note that the relations (17) and (18) are true whatever the axes, because 
using (8) and (9) , they imply Kjf, = Hj^. If we want that the amplitude 
remains unchanged, it is necessary that 

Hu = cos 9 

Here, the "amplitude" is the coordinate along the axis of polarization, and 
it is a quantity which can be complex. 

More generally, any beam is defined by two polarized waves, the first 
one in the direction Ox, the second one in the direction Oy. Each of them 
is submitted to a rotation of angle 9. The amplitude of each of them is the 
result of a LIF with complex gain H\\j cos 9 (which is a quantity invariant 
in a change of basis). Then we can speak about a rotation whatever the 
kind of wave, polarized or not. 
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Also we see that an unpolarized wave will be changed by the bi-filter in 
an unpolarized wave such that (when 9 7^ \ mod ir) 



\H 



11 



cos 2 9 



But more general conditions allow to retain the unpolarization. 



3.7 Depolarization 

The medium in free space or in optical fibers is a source of depolarization. 
Bi-filters take into account this situation. If we start from a polarized wave 
along Ox, we have by definition 



Kit) 

E Z y(t) 



7*21 



(<) 
(0 



Except when [H11/H21] (w) is a real constant, the wave is no longer polarized 
(assuming H\\ and H21 different of 0). With p x , p y , p z xy as defined in section 
3.2, we define the constants a, b, c, d by 



c 
d 



I- a 
-b = 



-e lv p xy /p' 



P Z xPy 



assuming that p' 7^ and where v is any real number. This leads to the 
Stokes decomposition [19] 

f E z x (t) = A(t) + B(t) 
< A(t) = aE* x (t) + bE*(t) 
( B{t) = cEl{t) + dE z y {t) 

with a polarized part (A, 0) and an unpolarized part (in the weak sense) 
(B, EfJ) , which verify the conditions 



E[B(t)E**(t)]=0, E \B(t)\ 



E 



\E z (t)\ 



To perform this decomposition with a given power for each part, it suffices 
to choose filters with complex gains Hn,H2i such that 



Py = 

p z x = 

Pxy 



00 

00 
00 

00 



I tt |2 n 
I-H21I s x 

\Hn\ 2 s x 



(lo) dio 
(w) doo 



The condition p' 7^ can always be verified using well-chosen filters. But 
other Stokes decompositions can be done, changing the definition or the 
basis [8]. 
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3.8 A particular class of beams 



The main problem is to know if the model is well-fitted to physical situa- 
tions. A particular model is treated in appendix 5. A quasi-monochromatic 
light can be polarized, or unpolarized or partially polarized. In the first 
case, the electric field has a constant direction and in the second case, the 
electric field takes any direction with equal probability (independently with 
its amplitude). They are extreme situations and the intermediary situation 
is for a partially polarized light where the field direction is a random process 
with a one-dimensional probability law which is not degenerate (such as po- 
larized light) neither uniformly distributed on (0, 2ir) (such as unpolarized 
light). A large class of light spectra verifies (see appendix 5) 

4H = a|/M| 2 

= (i-«)l/ Ml 2 (19) 
4» = /?|/H| 2 



where < a < 1, and \(3\ < y/a (1 — a). The values = \J a (1 — a) 
are for polarized light in a direction defined by a, and a = \,P = for 
unpolarized light. A quasi- monochromatic wave corresponds to / (oj) = 
outside a short interval centered at some frequency uo/2ir. (7) shows that 
the bi-filter does not change the membership to the class if the parameters 
Hjj, are constant on the spectral support of the wave. For instance, the wave 
(19) with parameters (a, P, /) is transformed in the wave with parameters 
(a', ft, g) by the bi-filter 

H22 = kH\2 



«n=*H„ <20) 

where the new parameters are (k, a, P are scalar, / and the Hjk are functions 
of oj) 



l + \k\ 2 J , p' = k* (l + \k\ 2 ^ 
\g\ 2 = \f\ 2 (l + |£f) ((1 - a) \H 12 \ 2 + a \H U \ 2 + 2pTZ [H 12 H^ 



4 About the Beer-Lambert law 
4.1 The Beer-Lambert law for niters 

In acoustics or ultrasonics, a beam is represented by a (real or complex) 
scalar quantity U z (t) parametrized by the distance z at a transmitter. Very 
often, the propagation medium is modelled by a LIF. If the medium can 
be split in "independent" pieces, the complex gain of the equivalent filter 
verifies the functional equation 

H Z1Z2 (w) = H Z1U (w) H UZ2 (to) (21) 
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where H ZlZ2 (oj) is the complex gain of the piece of medium in the interval 
(zi,z%) . (21) expresses that the pieces (zi,u) and (u, z-i) have the behavior 
of filters in series. Regular solutions of (21) on M + are in the form [6] 

But if we suppose that 

H uv (oj) = H v - U (oj) 

(this means that only the thickness of the pieces appears) the last equality 
is reduced to 

H z (oj) = e- za{u>) (22) 

where a (ui) depends only on the medium properties. For example, the 
model a (oj) = —aoj 2 (a > 0) is admitted in limited frequency bands for 
atmosphere or water acoustic propagation, but other functions can be used 
[2], [13], [16], [8]. For optical propagation in free space, the function a(oj) 
is very complicated due to deep absorption holes. (22) is the Beer-Lambert 
law used in many domains of science. Finally, the spectral density s z (oj) 
and the power P z of the process U 2 verify 

s z ( w ) = e -2*R[a( W )] fl ^ ^ P z = e -2*»[a(w)] a ^ doj . 



4.2 The Beer-Lambert law for bi-filters 

1) We ask the question to know if a generalization of the Beer-Lambert 
law can be done when bi-filters are used. Of course, we assume that the 
successive pieces of the medium do not interact and that their properties 
are independent of the beams which cross them. Each piece in the interval 
(u, v) is represented by four LIF of complex gains HJj! with respect to the 
basis Oxy. We assume some geometric coherence of the medium so that 

Hf k = H]-\ (23) 

Figure 2 gives the equivalent circuit for two successive layers. To alleviate 
the formulae, we generally omit the variable oj. Consequently to the (strong) 
condition (23), and using the scheme of figure 2, the problem can be trans- 
lated in a set of differential linear equations. We assume the existence of the 
derivatives h z k = -§^Hj k for z = 0. We prove in appendix 4 (case 1) that, 
when Ai, A2 are distinct different of roots of 

A 2 - {h° n + h° 22 ) A + h° n h° 22 - h° 12 h° 21 = (24) 

the only one solution of the problem is 



ttz _ fe ii~ A 2 ^zAi I A i- fe n „z\?, 

H z 2l = ^ 2 {e^-e z ^) 

22 Ai — A2 Ai — A2 



(25) 
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Figure 2: 
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These equations are for bi-filters the version of the Beer-Lambert law for 
filters where the medium is represented by the family of complex gains 
H z (ijj) = e ~ 2a ( w ). Also equations (25) define the set of admissible repre- 
sentations of a continuous medium with bi-filters in series. Firstly we study 
the case where the derivatives h z - k = §^H z k verify the conditions 

h°n = h° 22 , h° l2 = -h° 2l . (26) 

2) Let assume that h\ 2 does not depend on the frequency lu/2tt. We find 
that the complex gains H z k verify 

( TTZ TTZ TTZ TTZ 

I M ll - H 22i H 12 - ~ H 21 (2 7) 

\ Hfa = e zh °i sin zh1 2 , H\ x = e zh ^ cos zh\ 2 . K ' 

If h\ 2 is a real quantity, we recognize the formulas for a rotation of the 
polarized beam (sections 3.6 and 6.2). The angle of rotation is equal to 
—zh® 2 (proportional to z), and the amplitude A 2 is the result of a LIF 
with input A and complex gain e 2 ' 1 ' 1 ' 1 (^i can depend on uj). Then the 
Beer-Lambert law (22) is true for the amplitude when the relations (26) are 
fulfilled for real h\ 2 . Also is not real, because the imaginary part holds 
the propagation time between the points O and z (which cannot cancel and 
which depends on the frequency ui/2tt in case of dispersion). The real part 
measures the weakening of the wave between both points. 

To summarize, we can choose the parameters of the bi-filter to model a 
beam with a given polarization at 0, which is rotated by any angle, and with 
any weakening and with any propagation time. The (generally complex) 
amplitude obeys the usual Beer-Lambert law. The angle is ruled by h\ 2 , the 
weakening by TZ [h®-,] > the phase by I [^?i] , and the three parameters are 
proportional to the distance z. 

3) Now we consider any beam with its components E® (t) , E° y (t) . When 
(26) is verified for the h® k , these equalities are still true for the H z k , and 
these quantities are invariant by rotation {K z k = H z k using (9) and (10)). 
Then both components have the same behavior when crossing the medium, 
i.e. are rotated, delayed and weakened by same quantities. In the basis 
Ox'y' such that —zh!{ 2 = (Ox, Ox') , we have 

E% {t) = T [Eg] (t) , £ 2 , (t) = T [Eg] (t) 

where F (cj) is the complex gain of the LIF T . Consequently the power at z 
will be given by (from (8) and the Wiener-Lee relations, see appendix 1) 

P z = f°° e^KH] [si + s° y ] ( W ) (28) 

J — oo 

to be compared with the result in the standard case (see the end of section 
4.1). Then the power P z at z is independent of the rotation defined by the 
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parameter h\ 2 . Equivalently, when studying the power, the beam state of 
polarization is irrelevant, and the medium can be viewed as a simple LIF 
where the Beer-Lambert law is available with the parameter /t°i • 

However, if we leave off the hypothesis of independence of h\ 2 with the 
frequency, the interpretation of the bi-filter as a rotator is not true because 
the angle of rotation is now a quantity proper to each frequency. 

4) When h\ 2 = pa 1 ^ is no longer real (but still independent of the fre- 
quency), (27) is true but with complex trigonometric functions which are 
expanded in 



H 12 



a cos (zp cos (ji) + b cos (zp cos (p — § ) 
: a sin (zp cos <j>) + 6 sin (zp cos </> — §) 
a = e zh \\ cosh (zp sin cp) 
b = —ie zh it sinh (zp sin cp) . 



This means that a polarized wave is split in two parts, the first one rotated 
by the angle —zp coscp, with amplitude coming from a LIF A with complex 
gain a, the second one rotated by (—zpcoscf) + ^) and with the complex gain 
b of a LIF B. The power P z is given by (the calculus is performed in Ox'y' 
with —zpcoscf) =(Ox, Ox')) 

P Z =E{[\A [E°] - B [Eg] | 2 + \B [E°] + A [Ej] | 2 ] (t)} 

Computations can be performed (to simplify we give results for an unpolar- 
ized wave s® y = 0) 



P, 



e 2zn[h\,} ^ + s o^ CQsh ( 2p zsin0)j (w) dw. 



We see that (28) is no longer true. 

Assume that the wave is polarized at O. Two reasons lead to a depolar- 
ization of the wave. The first one is the dependency of the parameter h® 2 
on frequency. Then the angle of rotation of the beam is different following 
the frequency. The second one happens when /t°2 i s n °t rea l> which creates 
a secondary wave with amplitude b and orthogonal to the main wave (of 
amplitude a). Because the function tanh is increasing on R, we have the 
same property for the quotient \b/a\ (but which cannot reach 1). 

5) We go back to the general formulas (25) . They show that the beam 
can be split in two parts, the first one containing the terms with coefficient 
e zAi an( j foe second part with e zX2 . Generally Ai and A2 are not conjugate 
complex numbers because the equation (24) can have complex coefficients. 
Assume that the input is the pure wave e tuJot . If 



Ai = a + ib, A2 = c + id 
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the wave at z will be the sum of two pure waves which cross the medium 
with celerities —coo/b for the first one and —ujo/d for the second one and 
with attenuation ruled by a and c. Celerities are different when b ^ d as 
in a birefringent medium. Obviously, each wave obeys the Beer-Lambert 
law in its simplest form but not the sum, except when the hjk verify some 
conditions. Finally, equations (25) are given when the roots of (24) are 
distinct and different of 0. Other situations are described in the appendix 
4. 

5 Conclusion 

A bi-filter is defined by a circuit of four linear invariant filters. It is a par- 
ticular case 2x2 of the MIMO circuits (for multiple inputs-multiple outputs) 
which are used for instance in communications between systems of antennas 
[14j and in sampling to improve the reconstruction of signals [IT]. It gen- 
eralizes the well-known notion of "scattering matrix" in radar processing. 
In this paper, we address the problem of modelling a medium crossed by 
a beam with two components and with some degree of polarization. We 
assume that the beams are random processes with properties of stationarity 
and spectra with any bandwith. We prove that bi-filters explain elementary 
operations on electromagnetic waves and we establish a generalization of 
the Beer-Lambert law which justifies the model in continuous media. Ob- 
viously, other situations can be highlighted. The parameters of the bi-filter 
are defined by the physical properties of the medium. Perhaps theoretical 
considerations about the crossed material could allow the determination of 
these parameters but I feel that it is a difficult task. In the field of ultrason- 
ics, for instance, the attenuation and the celerity of waves in some medium 
(sea water, biological tissues...) are obtained by experiments and not from 
the mechanical, chemical. ...considerations. I believe that the same applies 
for electromagnetic waves. Characteristics of a material are measured by 
studying a set of waves with different frequency and polarization. These ex- 
periments are able to estimate the four complex gains which define a bi-filter 
and they can help to give a fair representation of the medium. 

6 Appendices 

6.1 Appendix 1: Bi-filter 

1) Firstly, we summarize the usual theory of linear invariant filtering (LIF) 
of stationary processes [1], [12], @j. If X = {X (t) ,t € M} is characterized 
by its spectral density sx (<^>) , it is possible to define an isometry Ix between 
the Hilbert spaces Hj and K s where 

-Hx is the set of linear combinations of the random variables X (t) 
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(completed by the closure of this set) when the scalar product (., .) H and 
the associated distance d are used: 



(X(u),X(v)) H = V[X 
d 2 (A,B) = E \\A-B\ 2 



(u)X* (v)\ 



-K s is the set of the / (w) (from R to C) such as \f (u) \ 2 sx {oo) duo < 
oo, with the scalar product (., .) K and distances defined by 



</,<7>K = I 

d 2 (f,g) = L 



oo 
oo 



fg*s x ] (uo)doo 



\f - g\ s x (oo)doo 
The isometry I x is defined by the relation 

X(t) <— ^ e iuit 



(29) 



The isometry maintains the scalar product and the distances of correspond- 
ing elements of the spaces. Consequently, it allows to perform calculations 
in the space K s rather than in Hj, using Fourier analysis and geometry of 
Hilbert spaces. 

If T and Q are two LIF with complex gains F (oo) and G (uo) , input X, 
outputs U= T [X] and \=Q [X] , we have 



( U(t)<—> Ix F (oo) e iut , V (t) <—> Ix G (lo) e 



iuit 



E [U (t) V*(t-T 
Su (uj)= \F\ 2 
{ s uv (lo) = [FG*sx}(oj) 



= J^ o e^r [ FG * Sx ]^ )duj 
(oo), s v (oo) = \G\ 2 s x (oo) 



(30) 



where sjj (oo) , sy (oo) , suv (<*0 are spectral and cross-spectral densities. 

Finally, if W=W [U] , and if H (oo) is the complex gain of the LIF H, we 
have 

Wi t) ^ Ix [HF] (,)e iul 

which is the relation for filters in series. 

2) In section 2, the processes Xi and Yi are stationary and stationary 
correlated. If we look at processes Y^ and Y" defined by 



Y X (t) = Y( (t) + Y" (t) 

y[ if) <— ^ 



f21£i 

S X , 



(oo)e 



iuit 



(31) 



Xi and Y 1 are the input and the output of a LIF filter of complex gain 
SYiXi/s Xl - Using (30) we obtain, whatever t, r £ R 



E[X, (t)Y['* (t-r)] =0. 
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This means that Y{ (t) is the orthogonal projection of Y\ (t) on Hxj and 
that Y" (t) is orthogonal to H Xl , and 



Y{(t)€H Xl Y{'(t)±U Xl 




\SYiXi 



I ,2 



where Sy^ , s XiY{ > S Y" ■ ■ ■ are spectral and cross-spectral densities, 
quently, (4) can be split in two orthogonal systems S' and S" 



(32) 



Chouse- 



s' 

s" 



(33) 



X> 2 (t) = Hu [Xi] (t) + H12 [Y[] (t) 
Y 2 ' (t) = H21 [Xi] (t) + "H 22 [Y'J (t) 

X 2 ' (t) = H 12 [Y'/] (t) 
y 2 " (t) = % 22 [Y'{] (t) 

Y l (t) = Y> it) + Y>> (t) 
Y 2 {t) = Y> {t) + Y 2 " (t) 
X 2 (t) = X 2 (t) + X%(t) 

By construction, the sets S' = (X 2 , Y 2 ) and S" = (X 2 , Y 2 ) are uncorrelated. 
Each equation of S' is equivalent to a circuit composed by three filters. X 2 

is the output of the LIF with complex gain Hu + S ^ lXl H\ 2 with input Xi. 

Sx i 

The filter of complex gain H 2 \ H — YlXl H 22 is for Y 2 . 



3) From (30) , (32) , (33) we deduce the spectral characteristics of (Yi, Y2) 
(we omit the variable uf) 



sx 2 Y 2 = (^ x fH 12 + H u ) {^H 22 + H 21 )*s Xl + H 12 H* 2 



Sy 2 



lfjV£il 



sx 2 = 

S Y2 = 
s XrX 2 



H\ 2 + Hu 



s Xl + \H 12 \ 



^H 22 + H 21 



SY 1 ~ 



s Xl + \H 22 
^H 12 + Hu)*s Xl 

SX^ = (^^^22 + H 21 ) 

s Yl x 2 = ^ (^^12 + Hu)*s Xl + H 
SY lY2 = S -^(^H 22 + H 2i y SXl+ H< 2 



\ S YlXl\ 
S X, 



s Yl 



(34) 



6.2 Appendix 2: Rotator 

1) We study bi-filters which have polarized beams as input and output. We 
look for the bi-filters which increase the angle of polarization by a given 
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quantity 9, independently of the absolute values of the angle at the input 
and independently of the power spectra. 

Assume that the beam is polarized at and z with angles ip and ip' = 
ip + 9 (with respect to Oxy). It is equivalent to have (section 3-1) 

A z (t) cos ip' = [7*n cos ip + H12 sin ip] [A ] (t) 
A z (t) sinV' = [H21 cos V + ^2 2 sinV] [A ] (t) 

Taking ip = and ip = ir/2 leads to 

H 21 [A ] =Hn [A ] tan# 
^22 [A ] = -H12 [A ] cote 

equalities which are true for any A . Then the equalities about the complex 
gains become 

H21 cos 9 = Hn sin 9 
H12 cos 9 = —H22 sin 9 



(35) 



We enter (35) in the first equality which becomes 

A z (t) cos (if; + 9) = [Hu cos ip - V-22 sin %j> tan 9} [A ] (t) 
A z (i) sin (if; + 9) = [Hu cos ^ tan 9 + ^22 sin ip] [A ] (t) 

and we deduce the equality 

. , n . Hu cos ib cos 9 — H22 sin -0 sin 9 

cot (ib + 9) = 

Hu cos V> sin + #22 sin ip cos 9 

which has to be true whatever ip. Obviously, it is possible if and only if 
#11 = H22- Using (35) , we conclude that a NSC for a bi-filter to induce a 
rotation of angle 9 is summarized by 

H21 = -H 12 = Hu tan 9 

H22 = Hu (36) 
I * (t) = [A°] (t) 

except 9 7^ \ mod7r. For this particular case, we have 

H\\ = #22 = 0, H\ 2 = —H21 
A z (t) = -H12 [A°] (t) . 

To summarize, bi-filters verifying (36) define transformations composed by 
a rotation of the direction of polarization, associated to a LIF for the am- 
plitude. The LIF has input A , output A z and complex gain Hu/ cos 9 (for 
9 + § modvr). 

Then a "pure rotation" which retains the amplitude corresponds to the 
bi-filter 

Hu = #22 = cos6> 
H 2 i = -H12 = sm.9 

which is not surprising. 



20 



6.3 Appendix 3: Unpolarized wave 

1) If s xy = whatever the system of coordinates, we have, when 9 
(Ox, Ox') 

s x'y' = { s y ~ s x) sin 9 cos 9 = 



s z x , = s z cos 2 9 + s z sin 2 9 (37) 

sf,i = sf. cos 2 9 + s z sin 2 9 

y y Jj 



whatever 9, and then 



q% G^ 

°x' ~ °y' ~ °x ~ s y 



The definition which we have taken for the unpolarized beam holds on the 
whole process E^ and not only on the two-dimensional random variable 
[E% it) , E z (t)) . It takes into account all [E% (t) , E z y (t')) . It is a strong 
difference. 

As an example, take E z (t) = % [E^] (t) , where % [..] is the Hilbert tran- 
form and E| is a real process. It is wellknown that this implies p z xy = 
0, p z = p z and these equalities remain true in any coordinates systems (the 
equations (37) are verified for variances and covariances). Also 

q Z 

a x ~ a y ~ a x ' ~ a y' 
S ly = S l'y< = -^^Sign 

where signu; = 1 for to > and — 1 for us < 0. Though p z y = in any 
system (s x (u) is even), the r.v. E z (t) and E z (t 1 ) are linked for different 
t, t' because the cross-spectrum is different of 0. 

2) A compensator corresponds to a bi-filter such that (in Oxy) 



H 



e~ iuje y 



9 X and 9 y are the delays applied to E x and E®. If we assume s xy = 0, the 
first formula of (7) implies that s z xy = 0. Then a compensator maintains the 
property of unpolarization when the strong definition is used. 

If we take E° = A + B, E° = A — B where B is the Hilbert tranform 

x y 

of A with real A, we obtain (sgnw = 1 for oj > and — 1 for uj < 0) 

s° x = s ° y = s z x = s z = 2s A , s° xy = -2is A sgn 

which implies p xy = because (w) is even. Then, the beam is unpolarized 
at O in the weak sense but not in the strong sense. However from (7) 

s z xy (u) = -2ie iu) ^y~ ex) s A {uj) sgnw 
Ply = 4 Jo°° S A ( w ) sinw (°y ~ °x) du. 

Obviously we do not generally have p xy = 0, which proves that a com- 
pensator does not maintain the unpolarization in the weak definition, ex- 
cept if the spectral support of sa (w) is small enough around some loq and 
cjq (9 y — 9 X ) close to a multiple of it. 
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Prom (7) and with the strong definition, the unpolarization is maintained 
from O to z if and only if 



H12H22 + HnH. 



ll- n 21 



0. 



6.4 Appendix 4: the Beer-Lambert law 

1) The beam state E^ = (E£,E*) at z is the result of the bi-filtering of 

by the Hjj; or the bi-filtering of E M ,u < z, by the HJ£ (of course they are 
function of oj but we can omit this variable). This leads to the equations 
(using elementary properties of circuits and looking at figure 2) 

tjOz JjOu TJUZ I tjQu TJUZ 

TjOz JjOu TJUZ I TjOu TJUZ 



'12 

TjOz 

H 21 



TJVU TJUZ i TJVU TJUZ 

ti \\ ti 2\ + ti 2\ ti 22 



which are simplified in (from (23) which translates the homogeneity of the 
medium) 



(38) 



-"11 


TJU TJZ- 


u 


1 TJU TJZ- 

+ -^21-^12 


U 


H 12 


TJU TJZ — 

- ti \2 ti \l 


u 


1 TJU TJZ — 

+ -^22-^12 


U 


TJZ 
H 21 


TJU TJZ — 

- -^11-^21 


U 


+ #21 #22 


U 


, H 22 


TJU TJZ — 


u 


1 TJU TJZ — 

t - n 22- n 22 


u 



For instance , we write the first equation under the form 

Tjz jja 1 

-"n tjz H n ~ 1 



Tjz+a 
H ll 



Tja 

, tjz n \2 

+ H 21 — . 



If we assume the existence of derivatives h 



jk 



m H jk we obtain 



h 
h 



12 

2 

21 



H i2 h n 



n 22 ~ 

with 



Hl 2 h 2 

h z — 
n jk ~ 



+ m 2 hf 2 
21 + -^21^22 



1 + H 22 h 22 



(39) 



_d_ TJZ 

dz n jk 



which includes the (realistic) conditions 

lim z ^o#Ti = lim^ #22 = 1 and lim^o-^12 = lim^ #2i = °- ( 40 ) 
The system can be written as the matricial equation 



h u = p H « p 








a 12 













hp 12 


n 21 





n 22 








hp 21 





^22 



(41) 
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Three cases can be highlighted, following the properties of P. We assume 
that Hj£ = because a wave in a unlimited medium is evanescent. This 
condition cancels constants which can appear in solutions of the system. 
2) We have three possibilities which are detailed below. 



Case 1 : 



(h°u ~ h° 22 ) 2 + 4h° 12 h° 21 + and h° u h° 22 + h° 12 h° 21 



H jk 



c jkl e XlZ + c jk2 e x * z 



where Ai,A2 are distinct eigenvalues of P (which have negative real parts 
for a passive medium). The conditions (40) lead to 





= d u e zXl 


+ (l-d u ) 


e zX 2 


( H 12 


= d 12 e zXi 


- d 12 e zX * 




H 21 


= d 21 e zX > 


- d 21 e zX * 




, H 22 


= d 22 e zXl 


+ (1 - d 22 ) 


e z\ 2 



(42) 



with, using (39) 



(in : 
d2i 



Ai— A2 : 



d\2 

d22 



h° 
Ai -A 2 



(43) 



A1-A2' ^ A1-A2 

where Ai and A2 are (distinct, different of and with negative real parts) 
roots of 

A 2 - {h° u + h° 22 ) A + h° u h° 22 - h° 12 h° 21 = (44) 

Conversely, (42) with (43) verify (38) . Moreover formulae (11) imply the 
invariance of Ai and A2 in any rotation. We obtain the set of K z - k fitted to 

Ox'y' by replacing h° jk by k° jk =^K^ k in (42) , (43). This case is developped 
in section 4.2. 



Case 2 : \ ^ " ^ + 4/l ^°i + and h° u h° 22 = h° l2 h° 2l 



H* = c 3k e Xz 



where A = h\ 1 +h [ 22 is the eigenvelue of P assumed different of 0. Conditions 
(40) imply 

en = c 2 2 = 1, C12 = C21 = 0. 

Obviously the usual Beer-Lambert law is verified. Each component is weak- 
ened and delayed through a quantity proportional to z. 

Case 3 : (h° n - h° 22 ) 2 + Ah° 12 h° 21 = . 

In the last case, P has only one eigenvalue (of order 4) wich corresponds to 
a proper subspace of dimension 2. Actually, we find same results as in the 
case 2. 
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Then we have shown that the system (38) has an unique solution most of 
the time depending of two parameters (for instance h^, h\ 2 )- However they 
may be functions of the frequency lo/2it. Consequently bi-filters are able 
to model propagation of electromagnetic beams through continuous (and 
stationary) media such as free space or optical fiber or coaxial cable. 

6.5 Appendix 5: a class of beams 

The simplest model of incoherent light is described by 

x (t) = J2j e^-^hit - tj) cos Qj 
Y (*) = Ej e^-^h (t - tj) sin 0,- 

where t={t n ,n € Z} is an homogeneous Poisson process with parameter A, 
the 9 n are random variables independent of t and between them, and h (t) 
is regular enough. Each term represents the emission by the particle j at 
the time tj in the direction Qj. Straigthforward calculations yield (assuming 
H (wo) = to suppress some continuous component) 

s x M = (w-w )| 2 E[cos 2 e] 
s Y (w) = ±\H (w- w ) | 2 E [sin 2 9] 
sxy (w) = \H (uj - w ) | 2 E [sin 6 cos 9] 
H ( w ) = JZo h («) e~ iuju du 

Taking 9 n uniformly distributed leads to (this means that the elementary 
emitters have no favourite polarization) 

f s x (u>) = s Y (w) = ^ \H (u - w ) | 2 

\ S X Y (w) = 

Then, the wave with components X, Y is unpolarized. 

When Pr[9 = 0] = 1, we have a polarized beam in the direction 9 : 

s x M = ^ |-ff (w - w ) | 2 cos 2 6* 
s r M = 1^ ( w ~ w o) | 2 sin 2 
«xr (w) = 4^7 l-ff - w ) | 2 sin cos 6 

Other laws for 9 give a large choice of situations in the form (for real a, j3) 

f s x (co) = a\f(co)\ 2 

< S y(w) = (l-a)|/(w)| 2 

{ Sxy (w) = /3|/(w)| 2 

where \/3\ < \J a (1 — a). Conversely, we can find a probability law for each 
value of (a, (3) . 
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